Air parcel random walk and droplet spectra broadening in clouds 
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We study the effect of turbulent flow on the droplet growth in a cloud during the condensation phase. Using 
the air parcel model, we describe analytically how the size distribution of droplets evolves at the different stages 
of parcel movement. We show that turbulent random walk superimposed on an accelerated ascent of the parcel 
makes the relative width of droplet distribution to grow initially as /'/^ and then decay as t^^l'^. 
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The growth mechanism of small water droplets in clouds 
is the condensation of water vapour Q]]. While the dynamics 
of droplet growth is well understood, no consistent statisti- 
cal theory exists yet for describing the evolution of droplet 
spectra in the cloud — see iH 01 and the references therein. 
The main difficulties in this task are that a consistent theory 
should describe the interaction of at least four fields in the 
cloud: vapour concentration, temperature, droplet distribution 
and velocity. For sufficiently small droplets, the concentration 
inhomogeneities due to droplet inertia are negligible 
In cloud cores, vapour concentration and temperature also can 
be considered smooth functions of position, so here we shall 
only consider their vertical dependence. It is the velocity 
field which has a turbulent nature and can be considered as 
a stochastic variable. The statistical description of random 
interacting fields belongs to the class of field-theoretical prob- 
lems which are generally unsolvable. In our case, neglecting 
the spatial fluctuations of concentrations and temperature, one 
can turn the problem into that of a random walk superimposed 
on a regular drift in the prescribed environment. Three limit- 
ing cases of that problem (with no drift, permanent-velocity 
drift and permanent-acceleration drift) are solved here analyt- 
ically. This might provide a useful opportunity for comparison 
with the data of observations and numerical simulations. We 
use the concept of an air parcel fT], which means that we con- 
sider a macroscopic volume of air, which is moving through 
the cloud with the droplets within it. Together with the deter- 
ministic motion and growth, the parcel participates in a turbu- 
lent random walk which can be modelled as brownian motion 
on a timescale larger than the velocity correlation time. We 
examine how the form of the droplet size distribution changes 
with time. It is shown that the dynamics of the form of the 
distribution (spectrum) strongly depends on the deterministic 
motion. The relative spectral width decreases while the parcel 
is moving upwards either with constant acceleration or with a 
constant velocity. When there is no vertical drift (the parcel 
comes into mechanical equilibrium with its environment and 
wanders around a fixed height), the relative width of the dis- 
tribution begins to grow. The form of the spectrum increases 
in width. We also show, that the distribution of droplets over 
sizes is asymmetric: power law at small sizes part is deter- 
mined by the activation of new nuclei while the large sizes 



part is exponential. 

Advection and diffusion of passive fields in turbulent fluid 
is usually described with the Kraichnan model (see e.g. f^). 
According to this model the velocity field can be considered 
a stochastic Markovian process. The stochastic description 
is valid when the fluid motion is chaotic. This is true in the 
case of fully developed turbulence which corresponds to large 
values of Reynolds number j^l- In real clouds the observed 
Reynolds number is of order Re ~ 10^ -f- 10^ 1 1], which is 
certainly enough to assume air velocity field stochastic. The 
validity of the zero correlation time assumption is proved by 
the following arguments: correlation times of velocity fluc- 
tuations do not exceed 10^ s, while the duration of the con- 
densational stage lasts lO-' s or more |1]. Thus a more re- 
alistic model, which takes into account the finite correlation 
times would give a small correction not exceeding 10%. In 
this paper we focus only on the vertical motion of a particle, 
because horizontal inhomogeneities have fluctuating nature, 
and will not lead to any new qualitative effect. We also ne- 
glect collisions and coalescence of droplets by considering the 
initial stage of growth when droplets are small and the only 
relevant droplet growth mechanism is condensation of water 
vapourl 1|. Using all these assumptions, one can write the set 
of equations describing the dynamics of a single droplet |1|,|2]: 



- = u{t) + v{t). (2) 

Here a is the droplet radius, A = is a variable, propor- 
tional to the droplet surface area, and s is the local super- 
saturation rate defined as s — n/no — 1, where n and mq are 
the local vapour concentration and the concentration of satu- 
rated vapour correspondingly. We have set the coefficient of 
propotionality to 1 , which can be always achieved by using the 
proper dimensional units. For any given droplet, the supersat- 
uration rate changes with time because the air parcel contain- 
ing this droplet is lifted up to the cloud and is being cooled 
along the way. The vertical velocity of a parcel (|2j consists of 
two components: the deterministic u{t), which corresponds 
to the macroscopic updraft flow and the random v{t) which 
is due to turbulent velocity fluctuations. We will consider the 
random velocity to be gaussian and short correlated 
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{v{tMt'))=Db{t-t'). 



(3) 



In the following we will analyze the dynamics of the droplet 
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FIG. 1: An example of the parcel deterministic velocity trajectory together with the relative and absolute distribution width for arbitraty chosen 
parameters to,uo,r. The three curves in each plot correspond to the different cases: parcel with constant acceleration (solid line), parcel with 
constant velocity (dashed line) and parcel with zero velocity (dotted line). 



distribution P{A,z,t). The quantity P{A,z,t)dAdz represents 
the fraction of the droplets which are located in the area 
(z, z + dz) and have the surface area in the interval (A , A + dA ) . 
As long as we do not take into account the interaction of the 
droplets, the dynamical equation governing the evolution of 
the spectrum is described by the Liouville equation that cor- 
responds to the single droplet dynamics J1I2> : 

{d,+s{z,t)dA + [u{t) + v{t)]d,}P{A,zJ)=0. (4) 

Applying the standard procedure of averaging over the white 
noise (|3} one obtains the Fokker-Planck equation: 

{d, + s{z,t)dA + u{t)d, - ^Ddl} P{A,z,t) = 0. (5) 

The second and third terms describe the advection of the 
droplets in phase space due to deterministic terms in equa- 
tions iU2\ . They are responsible for the motion of the distri- 
bution center. The third term comes from the stochastic part 
of the velocity and its main effect is the spectrum broadening. 
We assume that the initial droplet size is small so that we can 
set P{A,z,0) — 5(A)5(z). This equation can not be solved for 
general supersaturationfield.s(z,f). Howeverithas an analytic 
solution for a stationary linear profile s{z,t) — so+siz which 
we consider as a reasonable local approximation of the real 
supersaturation field. The expression for P{A,z,t) is rather 
bulky. However the size distribution P (A, f) = / (izP{A,z,t) is 
described by the simple expression 
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(6) 



HereAo(f) — JqAi' {so + siZp{t')) andzp(f|is a a deterministic 
parcel trajectory: Zp{t) — /odf'M(f'). The distribution width is 
5A2(f)=Di^fV3. Usually one is interested in the distribution 
of radius a; from (|6j we see that it will be substantially non- 
gaussian in our case as Pa{a,t) ~ 2aP{a^,t). Nevertheless, it 
will still have the bell-form and can be described by it's center 
ao{t) — \/Ao and width 
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One can see that in some situations the distribution in a can 
be narrowed, even if it is broadened for the variable A. The 



growth rate of the distribution width can be described by the 
variable ^ = 5A /Aq. The relative width of the distribution over 
radii is then 
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The r.h.s. of this relation is a monotonic function, so the rel- 
ative width of the distribution over a behaves qualitatively 
in the same way as the variable t,. Even the regular (non- 
turbulent) part of the parcel motion in a real cloud can be very 
complicated. At the beginning the parcel is being lifted up the 
cloud because of the convection mechanism. The parcel is ac- 
celerated until the velocity reaches some stationary value due 
to the drag force. During the motion the temperature of the 
parcel changes because of droplet growth and heat exchange. 
When the parcel temperature becomes equal to the ambient 
temperature, its velocity decreases to zero and it stays at this 
steady height. One can analyze such a movement by consid- 
ering three limiting cases: when the parcel is accelerated so 
u{t) = uo+at; when the parcel's velocity is constant: u{t) =u; 
and when the parcel has reached its steady height, and its ve- 
locity is zero, u{t) — 0. In the following sections we will an- 
alyze these three cases and show that they correspond to the 
different types of behaviour of the size distribution function. 

We start from the simplest case of zero velocity. As already 
mentioned, this case is important when analyzing the last part 
of the parcel's movement. On this part the distribution center 
has akeady moved from A = and its width is also finite. 
However if the steady stage lasts for a long time, the initial 
conditions are fo|-gotten, and one can treat this problem yith 
the initial conditions A (0) = and 5A(0) ~ 0. The solution of 
this problem was introduced in the previous section (|6jl, and 
one obtains: Zp{t) = 0, Ao(f) = sot, dA^{t) = Ds\t^/3. Thus 
for the relative width one finds: 



(9) 



Therefore the distribution is being broadened with a rate t, oc 
f '/^ at this stage. Turning to the radius distribution, Pa{a,t), 
one can say that its relative width is increasing with the same 
rate 5a/a°^ f '/^ for small times t < s^/s\D, but as ^ be- 
comes large enough according to (|8} its growth rate changes 
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to ha/a cx: ^'/^ o= f'/^. The radius distribution is broadened 
slower than that of the surface area. It is worth noting that in 
stratiform clouds where the deterministic parcel movement is 
small the zero velocity case is most important. We see,that 
one should observe strongly broadened distributions there. 

Next we consider the constant velocity case: Zp{t) = ut, 
Ao{t) =Ao{0)+S()t + {siu/2)t^. Again, at large enough times, 
the distribution forgets about the initial conditions (which 
were formed during the acceleration stage) and the last term 
in Ao(f) will become dominant, so Aq = {siu/2)t^, 5A^ ~ 
{Ds\/3)t^. Therefore the relative distribution width is slowly 
decreasing: ^ f while the center of the distribution and 
its absolute width are growing with the rates ao = \/Ao ^ t. 

Finally, we consider a uniformly accelerating parcel. Here 
one finds s{t) — s{zp{t)) — SQ + siuot+siat^. If initial droplets 
sizes are small, we can easily find the dynamics of the center 
of drops surface area distribution, Ao(f ) = ^of + {siUo/2)t^ + 
(iia/6)f^. The relative width behaves like 
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At the very small times the relative width is growing with the 
rate ^ 0= f as in the zero velocity case. At large times, when 
the parcel has accelerated enough, the relative width tends to 
zero as ^ 0= f^^/^. Therefore, there is a time when the rel- 
ative width has a maximum. The position of the maximum 
is determined by the values of uo,a. It can be estimated as 

fmax ^ min sq/ (siuq), yGo/Jsia) . The value of the maxi- 
mum is in this case ^max ^ \/Dt^si/so- It is interesting, 
that the absolute width 5a ~ ^/Ao{\/l 1) ~ \/Ao^ is al- 
ways increasing, and approaches a constant asymptotic value 
at large times. The particular value of this width is a compli- 
cated function of all the parameters. 

The figure[2shows the qualitative behaviour of the absolute 
and relative width for the three analyzed cases. The pictures 
should not be considered as a prediction for some particular 
cloud, but as an example, which shows how different velocity 
trajectory can affect the distribution parameters. 

So far we have analyzed the situation where small droplets 
are injected into the cloud only once, at the beginning. In real 
clouds the situation is certainly different. New nuclei are ac- 
tivated all the time, and their activation surely affects the dis- 
tribution form. Nevertheless, we will show that this effect of 
new nuclei activation does not change the large droplet sizes 
part of the distribution. The droplet size distribution can be 
written as follows: 



n(a,f)=/ dt'no{t')Paia,t~-t'), 
Jo 



(11) 



where no{t)dt is the number of nuclei activated in the interval 
{t,t + dt). Luckily the theory of nuclei activation is well de- 
veloped 1 1]. The simplest model of nuclei activation predicts, 
that the activation rate is proportional to some power of local 
supersaturation rate, no{t) °^ s^^{t), and usually k > 0. The 
absolute width 5fl(f) of the distribution Pa{a,t) is a growing 
function, therefore the front of the real distribution will be al- 
ways determined by the droplets which began to grow at f = 0. 
Nevertheless the new droplets determine the form of the left 
part. Assuming the distribution to be narrow enough, one can 
have the following approximation to the small droplet part of 
the spectrum: n(flo — Afl,f) ^ no(f*), where t* is determined 
by the condition that Pa{a,t*) should be centered at a — Aa. 
For example for the acceleration stage 5(f) = sq + siat^ /2, 
and the center position is given by Ao(f) = sqI + siat^ /6 °^ 
t^. Thus flo(f) °= t^^^- From the definition of t* we find 
f — f* cx (flo — Afl)^/^, so t* al^^ — {ao — Aaf-I^ oc a^^^^Aa 
for small Ao <C flo- Recalling that no(0 °= ^, we find that 
n(flo — Afl,f) 0= (io + cAfl^) . We can see that the left part 
will have a rather wide power law tail distinct from a narrow 
exponential right tail. This is also an important result, which 
might be observed in real clouds. 



To conclude, we have shown that the evolution of size dis- 
tribution of droplet spectra can be divided into three different 
stages. At each stage, the relative width of the distribution 
behaves according to different power laws: it grows with a 
rate ^ f at the zero velocity stage, and decays with the 
laws ^ f^'/^,^ 0= f^^/^ on the stages with parcel constant 
velocity, constant acceleration respectively. The dynamics of 
real air parcels is more complicated, and the equation ^ al- 
lows one to estimate the distribution width according to our 
model. However, from physical reasoning we see that real par- 
cel dynamics can be approximated by three consequent stages, 
where the relative width asymptotics were computed analyti- 
cally. 



We have also shown, that the left part of droplet size spectra 
is mainly determined by the process of nuclei activation dur- 
ing the initial stage of parcel movement. In this case one will 
observe a power-law tail in the small droplets part of spec- 
trum, in contrary to exponential tail for large droplet sizes, 
which is determined by the parcel's velocity fluctuations. 
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